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Abstract  f  (  ' 

l 

Research  on  energy  storage  and  transfer  in  molecular  hydrdgen  and  its  isotopes  lead  to  the  discovery 
of  mid-infrared  lasers  in  H,1  and  H^2  .  Specifically,  the  laser  lines  in  Hj1  are  at  3.84,  3.77  and  3.71 
um;  in  H^2,  they  are  at  4.717  4.60,  and  4.52  urn.  These  have  b|ei\  reported  by  us  recentlyd  ,2,3) .  Ue 
conclude  that  these  lines  in  both  isotopes  belong  to  the  a3Zg  —  c3Ru  transitions.  For  H^2  there  is  still 
some  doubt^ince  several  of  the  spectroscopic  coefficients  for  t^ie  H^2 ,  c3Hu  state  have  not  been  deter¬ 
mined.  In  this  report v  we  describeSAthe  conditions  under  which-e)|e  lasers  occur  and  show  how  this  leads  to 
a  plausible  explanation  for  their  existence.  We  derive  several  plausible  sets  of  value  for  the  spec¬ 
troscopic  coefficients  for  the  H,,2  c3Hu  stat?V\^^ 

Introduction 

\ 

Research  on  energy  storage  and  transfer  in  molecular  hydrogen  and  deuterium  led  to  the  discovery  of 
mid-infrared  lasers  in  these  gases.  Specifically,  laser  lines  at  3.84  pm,  3.77  um,  and  3.71  um  in  hydro¬ 
gen  and  at  4.71  um,  4.60  um  and  4.52  um  iij  deuterium  were  identified  and  reported ( 1 , 2 , 3) .  Our  conclusions 
are  that  these  lines  are  part  of  the  a3Eg  -  c3IIu  series  for  both  isotopes  of  hydrogen.  For  normal 
hydrogen  this  seems  to  be  certain.  For  deuterium  there  is  still  some  doubt  since  several  of  the  spec¬ 
troscopic  coefficients  for  the  deuterium  c3IIu  state  have  not  been  determined(4) .  In  this  report  we 
describe  the  conditions  under  which  the  lasers  occur  and  show  how  this  leads  to  a  plausible  explanation 
for  their  existence:  A  tentative  set  of  spectroscopic  coefficients  for  the  deuterium  c3Jlu  state  are 
eiven. 

Laser  Apparatus: 

The  lasers  were  discovered  during  a  search  for  weak  or  induced  dipole  lasers  in  molecular  hydrogen, 
deuterium,  and  deuterium  hvdridef5).  A  detaile.d  description  of  the  laser  apparatus  has  been  previously 
reportedfl).  The  laser  basically  was  a  glass  tube  6.3  m  long  and  2.5  cm  diameter  with  adjustable  internal 
nirrors  which  were  highly  reflective  in  the  laser  bands  and  nearly  transparent  for  visible  light.  The 
mirrors  were  dielectric  coated  and  "tuned"  to  the  desired  infrared  bands  with  output  mirrors  having  80%  to 
90%  reflectivity  and  the  back  mirrors  having  greater  than  99%  reflectivity.  Cavity  alignment  was  verified 
by  a  helium-neon  laser. 

The  lasers  were  excited  by  a  four-stage  Marx  bank  at  erected  voltages  of  225  KV  to  320  KV.  Current 
limiting  resistors  were  placed  in  series  with  the  discharge  to  prevent  discharge  collapse  and  oscillation. 
Gas  pressures  were  varied  from  about  *5  torr  to  10  torr,  and  two  gas  temperatures  were  used:  ambient  ("v 
300eK)  and  cryogenic  (n.  100°K).  Variations  in  voltage,  series  resistance,  pas  pressure,  temperature  and 
gas  composition  all  influenced  the  laser's  output  characteristics.  Figure  1  is  a  sketch  of  the  laser 
apparatus.  Figure  2  is  a  photograph  of  the  laser  in  the  laboratory.  Figure  3  is  a  picture  of  the  high 
voltage  end  of  the  laser  showing  the  current  limiting  resistor,  peaking  gap,  corona  ring  and  cryogenic 
jacket . 

Data/Data  Reduction: 

The  following  parameters  were  measured  directly  or  inferred  from  direct  measurements:  1)  initial 
voltage,  2)  voltage  vs.  tine,  3)  current  v£.  time,  4)  gas  composition,  initial  temperature  and 
pressure,  5)  total  laser  output  vs.  time,  and  6)  spectrally  resolved  laser  output  vs.  time.  Figure  4 
shows  the  voltage,  current  and  detector  signal  in  a  typical  experiment.  Figure  5  shows  the  spectrally 
resolved  components  of  the  laser  pulse.  In  many  cases,  as  shown  here,  an  identifiable  time  delay  (n>  .2  p 
sec.)  occurs  between  the  onset  of  the  strong  line  at  3.77  um  and  the  other  two  lines.  Additionally,  in 
order  to  facilitate  data  reduction,  a  data  field  of  voltage,  current  and  initial  gas  pressure  was  obtained 
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Figure  1.  Apparatus  for  the  hydrogen  and  deuterium  triplet  lasers. 

Figure  2.  Hydrogen  triplet  laser  in 
the  laboratory. 

1)  Laser  tube  (b.fe  meter  length  x  2.5 
cm  diameter)  on  a  wooden  optical 
bench. 

2)  Output  window  on  sealed  adjustable 
mount . 

3)  HeNe  alignment  laser. 

4)  Four-stage  Marx  bank. 

5)  Corona  ring  for  high  voltage  end 
(detail  in  Figure  3). 

6)  10:1  voltage  divider. 

7)  Storage  oscilloscope. 

8)  3/4  meter  grating  spectrometer,  and 

9)  Cryogenic  gold  doped  germanium 
infrared  detector. 


for  hydrogen  at  ambient  temperature.  At  high  currents  (i.e.,  high  voltage  and/or  low  series  resistance) 
the  discharge  would  become  an  arc  which  placed  an  upper  limit  on  current.  At  very  low  currents  the 
excitation  was  insufficient  to  produce  laser  action. 

The  instantaneous  state  of  the  electric  discharge  was  determined  from  the  measured  parameters  and  tht 
electron  drift  velocities  published  by  Englehardt  and  Phelps(6)  or  more  recently  calculated  drift  veloc¬ 
ities  from  Taylor(7)  based  on  the  Boltaman  transport  code  developed  by  Carleton  and  Megill(8).  Compu¬ 
tation  of  variables  proceeded  as  follows: 

a)  Laser  action  occurred  only  during  the  glow  discharge.  This  type  of  discharge  was  estimated  to 
fill  approximately  2  cm  out  of  the  2.5  cm  bore  of  the  tube.  Discharge  cross  section  was  taken 
to  be  'v  3  cm2. 

b)  The  voltage  measured  included  both  the  potential  drop  in  the  gas  and  the  electrode  potential 
drop.  In  this  work  the  electrode  potential  drops  were  assumed  to  be  small  so,  to  a  first 
approximation,  the  entired  measured  potential  was  assumed  to  be  across  the  gas  column. 
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Figure  4.  Typical  voltage,  current  and  detector 
outputs  (from  oscilloscope  traces) 


D  1  iu  TIME  t/i 

Figure  5.  Spectrally  and  temporally  resolved 

lines  of  the  hydrogen  triplet  laser. 


With  the  voltage,  discharge  path,  and  gas  density  known,  values  for  the  electric  field  strength 
(volts/cm)  and  for  the  electric  field  per  particle,  F./N  volt  cm  could  he  determined.  Then 
using  transport  data  from  Englehardt (6)  or  Taylor(7)  the  electron  drift  velocity  could  he 
obtained  from  their  published  calculation  for  a  given  E/N  in  hydrogen  gas.  From  the  same 
sources  the  characteristics  electron  energy  (electron  volts)  was  also  obtained. 

From  the  measured  current,  known  discharge  cross  section,  and  the  electrons  drift  velocity,  e 
concentration  of  conduction  electrons  was  calculated. 

The  calculated  electron  concentration  and  the  characteristic  electron  energy  can  be  used  to 
calculate  an  electron  energy  distribution.  The  simplest  thing  at  this  point  is  to  assume  a 
Maxvel 1-Boltzman  energy/velocity  distribution  for  the  electron. 

Finally,  the  cross  section  data  from  BSse(9)  for  the  electron  collision  excitations  of  hydrogen 
from  ground  to  the  first  triplet  state  can  be  used  with  the  available  free  electron  concentra¬ 
tions  in  the  desired  energy  ranges  to  calculate  the  rate  of  production  of  a  Eg  and  c  Hu  states. 
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TABLE  I 


Cross  Section: 

,  2 
3  cm 

Length: 

6.6  m 

Volume: 

2.1  l 

HYDROGEN  TRIPLET  LASER  OPERATING  CHAKACTL RISTICS 

Series  Limiting  Resistor:  100  .. 

Gas  Density:  6.5  x  10^/cm^  (hydrogen) 
Initial  Pressure:  2  torr  (T,  ambient) 


TIME 
(y  sec) 

a 

V 

(volts) 

I 

(AMPS) 

E 

v/cm 

U) * 

cm/sec 

Ek+ 

e  volts 

Y 

cm 

P 

3 

watts /cm 

|  TV  Lt 

joules/ cm 

0 

8.E4 

0. 

121. 

1.88E-15 

1.8E7 

3.8 

0 

0 

0.0 

0.5  j 

7.6E4 

360. 

115. 

1. 78E-15 

1.7E7 

3.6 

4.22E13 

1.3E4 

6.  5E-3 

i.o  ! 

7.3E4 

520. 

111. 

1.72E-15 

1.6E7 

3.3 

6.5E13 

1.8E4 

1.6E-2 

1.5 

7.0E4 

580. 

106. 

1.64E-15 

1.5E7 

3.2 

7.7E13 

2.0E4 

2.6E-2 

2.0 

6.5E4 

520. 

99. 

1.53E-15 

1.4E7 

3.1 

7.4E13 

1.6E4 

3.4E-2 

2.5 

6.0E4 

420. 

91. 

1.41E-15 

1.3E7 

3.0 

6.4E13 

1.2E4 

4.0E-2 

5.6E4 

300. 

85. 

1.31E-15 

1.2E7 

2.8 

5.0E13 

0.8E4 

4.4E-2 

LASER:  Start  0.9  0  sec 

Peak  1.5  o  sec  (detector  output, 

6  volts),  (9.9  Kj/f-atmos  input) 
End  2.4  u  sec 


Data  Ref:  Case  B,  6  August  1982 

#  E4  *  10^,  etc. 

*  Electron  drift  velocity  -  current  carriers 
+  Electron  characteristic  energy  (electron 

"temperature") 


9583805Crm2l0f  Yv, P°,t.C"tlal  Wells  for  tbese  stfltes  are  95936.1  cm'1  for  the  "a”  state  and 

or  !i  pract  cari  nuroo  eSt:he‘  Y*  Y  #Bd  "C”  Bt3te  curves  ««  almost  identical; 

tor  ail  practical  purposes  they  have  the  same  Franck-Condon  characteristics.  We  therefore 

assume  that  they  are  nearly  equally  populated.  Table  I  shows  represent.* Iv*  samples  of  he 

current  and  voltage  data  obtained  and  reduced  according  to  the  above  technique. 

Gain  and  Laser  Output: 

butionVthe  rotational 'states"^  YhyY?  'T*'  Y*"  fr°m  the  ortho'Para  Population  distri- 

tllrfor  laser  gain?  ^^erium  hydride  there  is  no  ortho-para  asymmet,.  and  therefore,  no  poten- 

csnu  elec??oni?1!tl?ted  ^  Popul®tionf  for  varl°u®  vibration-rotation  combinations  of  the  o’lg*  and  the 
electronic  states  assuming  thermal  equilibrium  within  the  rotational  state  and  equal  total  state 
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populations.  Temperature  of  100°K  and  300°K  were  assumed.  Ar  300'K  in  hydrogen  the  favored  rotational 
transition  is  from  .1=1  to  J  =  2,  that  is,  this  transition  has  the  highest  population  ratio.  This 
matches  conveniently  the  hvdrogen  line  assignments  based  on  measured  and  calculated  wavelengths.  For 
deuterium  the  case  is  not  so  clear.  Our  simple  rotational  equilibrium  modelling  shows  that  the  J  *■  2  to  J 
«  3  transitions  are  somewhat  more  favored  than  the  J  ■  0  to  .1  '  I  transition  at  100°K  and  higher.  In 
spite  of  this,  it  appears  that  the  J  *  0  to  J  =  1  transitions  are  the  ones  that  are  running  in  the  deu¬ 
terium  laser. 

Calculated  and  Measured  Wavelengths: 

The  standard  spectroscopic  formulations  and  coefficients  for  energy  levels(4)  were  used  to  calculate 
the  vibration-rotation  state  energies  within  the  a3Ig  state  and  the  c3Hu  state  in  molecular  hydrogen. 
Table  II  shows  the  results  of  these  calculations.  By  following  the  A  J  =  +1  and  Av  =>  -1  rules  for  rota¬ 
tion  and  vibration  transition,  and  looking  first  at  the  favored  J  transition,  i.e.,  J  »  1  to  J  *  2,  one 
readily  finds  a  coupled  set  of  transitions  which  have  wavelengths  matching  those  observed  in  the  labo¬ 
ratory.  These  are  the  transitions  between  the  energy  levels  shown  in  Table  II. 


TABLE  II 


IDENTIFICATION  OF  THE  HYDROGEN  TRIPLET  LASER  FROM  ENERGY  LEVELS  CALCULATED  FROM  PUBLISHED 
SPECTROSCOPIC  DATA 


T 

e 

U) 

e 

U)  X 

e  e 

a3Ig+ 

95936.1 

2664.83 

71.65 

c3Hu 

95838.5 

2466.89 

63.51 

v  =  0 

TERM  VALUES 
a3Hg+,  J  =  1 

c3riu,  J  =  2 
^  97238. 

tu  y 

e  c 
0.920 


Av ( cm- 1 ) 
2601. 


B 

34.216 


1 (microns) 
3.84 


.  101782. 


104481. 


Deuterium  Energy  Levels: 


Assigning  lines  is  more  difficult  in  the  case  of  deuterium.  The  problem  which  was  straightforward 
for  hydrogen  must  be  worked  differently  for  deuterium  since  most  of  the  c3Ilu  state  coefficients  have  not 
been  determined.  We  must  infer  a  set  of  spectroscopic  coefficients  which  placed  the  observed  lines  in 
plausible  transitions.  In  standard  spectroscopic  notation(4)  the  bottom  of  the  potential  well  (electronic 
"zero"  point)  is  defined  by  T  ,  the  vibration  harmonics  spacing  by  w  ,  and  the  anharmonic  correction  by 
0!  x  and  w  v  (and  sometimes  fiigher  terms);  rotational  energy  is  characterized  bv  coefficients  B  and  a  . 
T^e‘  energveleve]  term  value,  T(cm  *),  to  this  degree  of  approximation  is:  ^  e 


T  +  u  (v  +  U)  -  u  x  (v  +  *s) 2  +  w  y  (v  +  *s)3 
e  e  e  e  e  e 


(8  -  a  (v  +  ’-i))  J(J  +  1) 

e  e 


w>here  v  and  .T  are  the  vibration  and  rotation  quantum  numbers  respectively.  In  the  case  of  deuterium 
the  coefficients  for  the  asFg  state  are  known  but  for  the  c3Ilu  state  values  of  T  ,  w  ,  u  x  ,  w  y  ,  a 
and  B  are  not  given(4).  The  data  we  have  are  the  three  measured  laser  lines  andetheepubl ished|4|  lowlst 
lying  energy  level  for  the  c3IIu  state.  With  this  limited  data  we  make  assumptions  as  follows  to  derive 
the  coefficients: 

a)  We  assumes  values  for  Be  s  15.  and  a  %  .5,  values  which  we  estimate  to  be  within  ±20%  of  the 
true  value.  Errors  in  these  coefficients  are  least  significant  for  J  »  0,  where  we  expect  to 
find  the  transitions. 

b)  From  the  published  vQ0  ,  the  lowest  lying  level,  we  obtain  a  term  value  for  the  c3Ilu  (v  =  0,  .1 
■  0)  state  of  96729  cm-1  . 

c)  We  have  assumed  the  (v  +  ij)3  term,  w  y  ,  is  the  highest  correction  term  required  so  that  the 
derivative  d3T/d(v  +  '-})3  is  a  constant!  6  w  y  ,  and  all  higher  derivatives  are  zero. 

d)  We  assume  a  kinetics  relationship  identical  to  that  of  hydrogen,  i.e.,  the  three  lines  are  in  a 
vibrationally  coupled  cascade. 

The  problem  now  may  he  solved  by  solutions  of  simultaneous  linear  equations  in  which  the  values  for  T  ,  u 
u>  x  and  w  y  are  the  unknowns.  Another  method  is  to  solve  the  problem  by  inspection  and  trial  solutions? 
In  this  method  the  equations  for  T  and  its  first  three  derivatives  are  treated  as  difference  equations 
and  term  values  are  generated  by  assuming  three  adjacent  vibrational  transitions  in  the  J  ■  0  to  J  *  1 
series  and  in  the  J  “  2  to  J  *  3  series.  The  objective  is  to  derive  values  for  T  ,  u  ,  <d  x  ,  and  u>  y 
which  are  plausible,  i.e.,  Te(c)  should  be  about  100  cm  1  lower  than  T  (a),  id  for  "c"  ana  "a"  states 
should  be  reasonably  close  together,  <d  x  should  be  a  positive  number  In  the  range  of,  say,  30  to  60  cm  * 
and  ui  y  should  preferably  be  approximately  *1.  These  "requirements"  only  say  that  the  c3IIu  state  of 
deuterium  is  similar  to  its  neighboring  "a"  state  and  the  whole  structure  is  analogous  to  the  like 
hydrogen  states.  We  have  tried  a  number  of  combinations  the  most  likely  two  of  which  is  shown  in  Table 
IIIA:  Table  IIIB  shows  an  alternative  fit  which  is  almost  as  good.  The  best  fit  to  the  data  appears  to  be 
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for  the  vibrational  cascade  starting  with  v  =  6  and  .1  =  0  as  suggested  to  nr 

Herzhergf 1 1 ) .  In  this  case  we  determined  that  term  coefficients  for  the 
1 7f\Q  /.i  v  >  /.n  cc  - i  , ,  . .  _  t  it 


i'-.  !\ihrowski  .mr* 


Hu  state  are:  T  =  <<S^V 


!?!k’  ‘"eXe-\‘i0'88,/ud  V*  -,l-2S-  These  vaU>‘“!»  Plus  our  assumed  ones  for  h'  and  a  generate  lire  wave 
-e  Vilhin/  C*  *  °7ho*e  derived  experimentally.  Our  values  mult  he  considered  as  tent" 
nllu  fht  he^  are  br!ed  °r  f°Ur  data  polntK  and  show  only  that  our  data  ran  he  used  to  produce 

wou?d  h^r  ^  ’  Went ificat ion  of  additional  line  in  this  cascade,  particularly  of  longer  wavelengths 

III  f6Sh  a  qUeFtl°n  °f,Where  the  cascade  starts  and  other  lines  at  different  temperatures 

would  be  valuable  in  determining  J  dependence.  ' 
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TABLE  IIIA.  IDENTIFICATION  OF  OUR  BEST  FIT  TO  THE  DEUTERIUM  TRIPLET  LASER  DATA 


T 

e 

!  w  i 

1 _ -JS _ j 

Ui  X 

e  e 

u  y 
e  e 

B 

e 

U 

3*.  + 
a  4 

g 

95958.1 

1885.84 

35.96 

0.340 

17.109  1 

0.606 

if  3rr 
c  Hu 

95855. 

i  1769.0  ! 

40.88 

1.25 

15.* 

.5* 

TERM  VALUES 


a3Ig+,  J  =  0 


c  Hu,  j  =  i 


101626. 


Av(cm_*) 


X (microns) 


2121. 

<i .  71 

2172. 

4.60 

2210. 

4.52 

*  Estimate  only  -  no  temperature  variation  data  obtained. 

H  Our  proposed  values  based  on  transitions  suggested  by  I.  Dabrowski  and  C.  Herzberg 
(See  Ref.  3).  a 


TABLE  I I IB. 


AN  ALTERNATIVE  FIT  TO  THE  DEUTERIUM  TRIPLET  LASER  DATA 


T 

a; 

1  CJ  x 

1  w  y 

-  , 

a 

e 

e 

1  e  e 

e  e 

e  | 

e 

3-  + 

3  L 

g 


96958.1 


1885. 84 


35.96 


0.340 


17.109  0.606 


3r 

c  Ha 

95892. 

1682.4 

15.75  -0.83 

15* 

.5* 

TERM  VALUES 

n 

. 

3  + 

,  J  *  0 

c3FIu,  J  =1 

*  Av(cm  3) 

[  X (mi crons) 

V  =  2 

100014. 

2119. 

4.72 

3 

1021. 33,  “ 

101578. 

2169. 

4.61 

4 

103747. 

J»10 3094. 

2205. 

4 . 54 

5 


105299. 


*  Estimate  only  -  no  temperature  variation  data  obtained. 
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